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Abstract 

Spin-transfer is a typical spintronics effect that allows a ferromagnetic layer to be switched by 
spin-injection. Most of the experimental results about spin transfer (quasi-static hysteresis loops 
or AC resonance measurements) are described on the basis of the Landau-Lifshitz-Gilbert equation 
of the magnetization, in which additional current-dependent damping factors are added, and can 
be positive or negative. The origin of the damping can be investigated further by performing 
stochastic experiments, like one shot relaxation experiments under spin-injection in the activation 
regime of the magnetization. In this regime, the Neel-Brown activation law is observed which 
leads to the introduction of a current-dependent effective temperature. In order to justify the 
introduction of these counterintuitive parameters (effective temperature and negative damping), a 
detailed thermokinetic analysis of the different sub-systems involved is performed. We propose a 
thermokinetic description of the different forms of energy exchanged between the electric and the 
ferromagnetic sub-systems at a Normal/Ferromagnetic junction. 

The derivation of the Fokker-Planck equation in the framework of the thermokinetic theory allows 
the damping parameters to be defined from the entropy variation and refined with the Onsager 
reciprocity relations and symmetry properties of the magnetic system. The contribution of the 
spin-polarized current is introduced as an external source term in the conservation laws of the 
ferromagnetic layer. Due to the relaxation time separation, this contribution can be reduced to an 
effective damping. The flux of energy transferred between the ferromagnet and the spin-polarized 
current can be positive or negative, depending on spin accumulation configuration. The effective 
temperature is deduced in the activation (stationary) regime, providing that the relaxation time 
that couples the magnetization to the spin-polarized current is shorter than the relaxation to the 
lattice. 

PACS numbers: 72.25.Hg, 75.47.De, 75.40. Gb 
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In the context of spintronics, the electrical resistance of magnetic nanostructures are 
tuned with the magnetization states. Giant magnetoresistance (GMR), or anisotropic mag- 
netoresistance (AMR) allow the magnetization states of nano-layers to be measured with 
great precision. Such magnetoresistances are easily scalable reading processes and are used 
for magnetic sensors and random access memorie (MRAM) technology. The possibility 
of controlling the magnetic configuration of a magnetic nanostructure by injecting spins 
emerged only in recent studies, opening the way to a readily scalable writing process for 
MRAMs application. This approach is also extended to thermally assisted switching, in 
which the heat fluxes are also exploited in order to help the magnetization reversal. In 
order to control the magnetic configurations and their stabilities (for reading and writing 
processes), in such magnetic nanopillars, it is necessary to understand on one hand the pro- 
cesses responsible for the magnetization reversal (in the presence of a magnetic field and 
heat), and on the other hand, the processes governing spin-dependent electronic transport 
at normal/Ferromagnetic interfaces. Taken separately, both effects are rather well under- 
stood today. However, coupling the two processes leads one to consider a large variety of 
possible mechanisms, called spin-transfer, that involve an ensemble of non-equilibrium sub- 
systems in interaction, with different populations of electrons and different populations of 
spins. The present work tries to clarify this picture with a phenomenological analysis based 
on non-equilibrium thermodynamics of open systems. 

Magnetization reversal provoked by spin injection has been observed in magnetic nanos- 
Tuctures of various morphologies, from spin-valve multilayers ;o nanowires 

contacts jlll . 12, [l^, l^j], and different types of magnetic domain walls 
. In order to describe and interpret these observations, physicists 
where forced to add one or two current-dependent terms into the well-known dynamical 
equations that describe a ferromagnetic layer coupled to a heat bath (Fokker-Planck or cor- 
responding Landau-Lifshitz-Gilbert equations). However, the question remains open about 
the deterministic (e.g. spin-torque) or stochastic (e.g. irreversible) nature of the terms to 
be added. 

It has been observed that for a time window larger than the nanosecond time scale, and in 
the framework of one-shot measurements (i.e. non-averaged, or irreversible measurements), 
the magnetization reversal induced by spin-injection is an activated process, with two level 
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fluctuations 
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22L I23L l24j or simple irreversible jumps [22, 25]. In these experiments, 
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governed by stochastic fluctuations and noise, the observed effect is accounted for by a cur- 



rent dependent effective temperature in the Neel-Brown activation law [22j. In contrast, 
for quasi-static measurements (e.g. magnetoresistance measured as a function of the mag- 
netic field or current with DC systems or lock-in detection system) and for high frequency 
measurements, oscillations and resonances indicate, in the frequency domain, the manifes- 



tation of quasi-ballistic precession effects 



iim 
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271 |28| . In these last experiments, 



the stochastic nature of the signal is averaged out, and the behavior is described in terms 
of current dependent negative damping within a generalized Landau-Lifshitz-Gilbert (LLG) 
equation. This negative damping formulation is motivated by the pioneering works of Berger 
2t| and Slonczewski j^] about the deterministic spin transfer torque theory. However, the 
deterministic approach cannot directly account for the magnetic relaxation measurements 
performed in the activation regime (as discussed in Sec. IV-A below). The hypothesis of the 
Slonczewski's spin-torque (presented as a current-dependent deterministic term in the mi- 
croscopic Landau-Lifshitz-Gilbert equation) is not useful as such in the description proposed 
here, i.e. in the context of open systems. 

In order to justify the introduction of the counterintuitive phenomenological parameters 



(effective temperature and negative damping), a detailed analysis o 



is performed on the basis of thermokinetic theory 
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44j . The first step (first section below) is to identify the relevant sub-systems of 



interest ( pointing out the difference between the spin-accumulation due to the diffusion 
of spin-dependent conduction electrons at an interface, and the magnetization of a ferro- 
magnetic layer), the coupling between them, and the role of microscopic degree of freedom 
that will be reduced to the action of the environment. In section two, spin-injection and 
spin-dependent transport are described in the framework of the two spin-channel approxi- 
mation (a conduction channel that carries spin up and a conduction channel that carries spin 
down, defined by the conductivities). Giant magnetoresistance, spin-accumulation, and cor- 
responding entropy production, or heat transfer, are deduced. Beyond the two spin channel 
approximation, the analysis is extended to four channels with the introduction of two other 
electronic populations (typically s-like for conduction electrons, and d-like for the ferromag- 
netic order parameter) and the relaxation between them. In the same manner as spin-flip 
scattering coupled the spin-up and spin-down channels, this relaxation defines a dissipative 
coupling between the ferromagnet and the spin-dependent electric sub-systems. The third 
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section is devoted to the detailed description of the ferromagnetic order parameter coupled 
to a heat bath (without spin- injection). Both the rotational Fokker-Planck equation and the 
corresponding LLG equation are derived in the framework of the thermokinetic theory, i.e. 
with the help of the first two laws of thermodynamics and the Onsager reciprocity relations 
only. The coupling of the ferromagnetic order parameter to the heat bath is introduced 



via the chemical 
temperature 



31 



otential with a typical Maxwell-Boltzmann diffusion term including the 
. l38j|. The Neel-Brown law is deduced in the activation regime. 
The last section is devoted to the ferromagnetic Brownian motion activated by spin- 
injection. The contribution of the spin-polarized current is introduced by the s — d like 
relaxation, as a source term into the conservation laws of the magnetization. Explicitly, it is 
shown that if n is the density of magnetic moments oriented in a given direction 9, $ of the 
unit sphere, and Jm is the corresponding flux of magnetic moments (this flux is not a dis- 
placement in the usual space), the conservation of n writes: dn/dt = —<HvJm+ f N _ F ^(z)dz, 
where the divergence is defined on the sphere and \1/ is the relaxation rate, integrated through 
the Normal- Ferromagnetic interfaces. This equation defines the irreversible spin-transfer oc- 
curring in the ferromagnetic layer, taken as an open system. The relaxation rate is related 
to the spin-accumulation Afj, through an Onsager transport coefficient L, \1/ = LAfi (where 
A/x is proportional to the current). L is linked to the relaxation times through the charge 
conservation laws (or electric screening properties). 

Due to the large relaxation time separation, the contribution of the source term can be 
reduced to the effect of an environment that is responsible for an effective damping and 
effective fluctuations (or effective temperature). The energy transferred between ferromag- 
netic layer and the sub-system defined by the spin-accumulation conduction electrons can 
be positive or negative, depending on the sign of the spin accumulation at the different in- 
terfaces. The effective temperature is deduced in the activation (stationary) regime, because 
the relaxation time that couples the magnetization to the spin-polarized current short cuts 
the relaxation to the lattice. 
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I. THERMOKINETIC APPROACH 



A. Interacting sub-systems 



35 



The general scheme of the thermokinetic approach is described in the references |32l 134 . 
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38| . The method consists in defining the state of the system with a set of the relevant 
extensive variables, say {s,Xi}, where is, e.g. the densities of particles in the sub-system 
i, or equivalently, the density of component i of a multicomponent fluid, and s is the total 
entropy density. The conservation equations should then be written, and the two laws of 
thermodynamics applied. The conservation equation for the component i writes: 

dn ' — * 

= -div(Ji) + Y.ju.^j (1) 

The divergence of the current Jj describes the conservative part of the process, and 
the term Sfrj is a source term that describes the relaxation of Uij components i into the 
component j iv^ < 0), or inversely [y^ > 0) 45j] . It is proportional to the inverse of the 
relaxation time tyj oc r _1 (see Appendix A). Physically, the term tyj describes the relaxation 
process that changes the internal degree of freedom (e.g. spins, electric charges, internal 
configuration). In terms of chemical reactions, tyj is the velocity of the reaction, i.e. the 
generalized flux thermodynamically conjugated to the chemical affinity Aj (defined below). 
The summation over all sub-systems, or all components of the fluid is that of a conserved 
variable: = — div(EiJi). The same holds, of course, for the energy E: ^ = —cIw(Je), 

where Je is the flux of energy. In contrast, the entropy production of the total system is 
not conservative in general, due to the irreversible processes (in other terms, information 
is lost). The equation for the entropy production of the whole system takes the canonical 
form || = —div(J s ) + X, where J s is the flux of entropy, and X is the internal entropy 
production, or irreversibility, which is a consequence of the second law of thermodynamics: 
X > (assuming T > 0). According to the first law of thermodynamics, the energy E, is a 
state function that is also scalar, extensive and conserved, so that 

dEjs^x,}) _dEds BE dxj 

dt ~ ds dt l dxi dt [ } 

where dE/ds = T is the temperature, dE/dxi = Fi is the generalized force associated 
with the flux dxi/dt. In the following we will deal exclusively with the chemical potentials 



6 



fj,i = dni/dE, unless specified otherwise (i.e. there is no need to introduce other extensive 
variables). The following Gibbs relation is obtained as a direct consequence of the first law: 
T|| = ^3 — Sj^-^Xj. After having inserted the conservation equations, Eq. (JTJ), into Eq. 
(J2J), the following form is obtained Q : 

ds -> -> 

T- = -div(J E ) + Z m diviJi) - E im (3) 

Using the development div{ni Ji) = ^ div(Ji) + Ji.grad(fii), Eq. (jSJ) can be re-written in 
the canonical form: 

?l = -div(J.)+I (4) 



where: 



1 = J E - grad C^J - EiJi ■ grad (J^j - //; v^j 



where the last term on the right hand side defines the dissipative coupling between the 
sub-systems. As will be shown in the last section, this term is responsible for the irreversible 
spin-transfer effect described in this work. What is unusual in dealing with the second law, is 
to manipulate an inequality instead of an equality, and consequently to deal with sufficient 
conditions instead of equivalences. Here, the condition X > leads to a positive matrix 
{Cij}ij of Onsager-Casimir transport coefficients that are state functions of the variables 
{s,Xi}, in order to build a positive quadratic form. The condition is fulfilled if the flux Ji 
and the relaxation velocity have the form 



— * — * 

Ji = — T>j£ij grad(jjj) 



— ^- J jJ J ij Aj 



(6) 



where 



Aj = —^Vik (7) 



is the chemical affinity of the corresponding reaction j (and we have Aj = —dE/d^j) |4j 
Furthermore, due to the time reversal symmetry of the microscopic equations, the transport 
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coefficients follow the Onsager-Casimir reciprocity relations [481 ] . The cross-coefficients that 
couple the flux J to the relaxation process \I/ are assumed to be zero, because, according to 
the Curie principle, only processes of identical tensorial nature are coupled. Inserting Eq. 
© into the continuity equation Eq. (Q), we obtain an equation of the time variation of the 
density drii/dt in terms of derivatives of the chemical potentials /i, : 

= ^j^-ij V fij + ^jk^ijLik A k (8) 

It is then sufficient to know the form of the chemical potential as a function of the density 
(for pur fluids : //(rij) = fio + kT ln(rii/N)) in order to derive the corresponding differential 
equation, or Fokker-Planck equation, with diffusion and relaxation terms (see sections II, 
III and VI below). 

What we gain in performing this analysis is to identify clearly the conservative and dissi- 
pative flux (through the internal entropy production), and to be able to define a dissipative 
process that couples the sub-systems beyond the usual deterministic coupling (electric field, 
magnetic field, pressure, etc. . . ). This dissipative coupling appears with an additional trans- 
port coefficient L, defined univocally via the transport equations. In the case studied below, 
the matrix C is composed by theficonductivities cxj associated to each channel (i.e. associated 
to a given electronic population), the thermal conductivity, or the corresponding Seebeck 



(thermoelectric power) and Peltier coefficients |43|, |49|, |50 1 and the ferromagnetic transport 
coefficients: gyromagnetic ratio V and the Gilbert damping coefficient 77. Beyond, the flux 
of entropy or heat allows the spin transfer to be understood in an open system in terms of 
relaxation with a supplementary Onsager coefficient L. As shown in the last section, this 
term is responsible for an effective temperature T e ff and effective (negative) damping a e ff. 

B. The model 

The model is based on the hypothesis that the ferromagnetic order parameter M is 
well differentiated from the sub-system composed by spin-polarized conduction electrons, 
although both systems exchange charges, spins, and heat through a relaxation mechanism 
that will be described in terms of internal variables As shown above, the relax- 

ation of an internal variable (or internal degree of freedom) defines a transport coefficient 
L s d related to the corresponding relaxation time r s d (L s d oc r^ 1 , see appendix A for the 
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relation to the relaxation time). 

We hence start with the two sub-systems: the ferromagnet described by the magnetization 
M and the two conducting spin-channel system of the conduction electrons. Both sub- 
systems are dynamically coupled through the relaxation time T s d- This relaxation is qualified 
as interband relaxation, to be opposed to the intraband spin-flip relaxation r s f introduced in 
the usual two spin-channel approximation. The conducting channels are usually described 
by the density n-f of conduction electrons with spin up and the density nj, of conduction 
electrons with spin down. The intraband coupling (accounted for by L s f or r s f) is responsible 
for the spin-accumulation mechanism at stationary regime. For convenience, we redefine the 
two channels with the density of spin-polarized electrons An = — nj_ ("spin conduction 
channel") and the total density of electrons no — n-f + n^. 

Furthermore, the conduction channels are contacted to a power supply (current generator 
here). Strictly speaking, the magnetic system is also contacted to the power supply, e.g. 
through the electron of d character j^jj. The conduction electrons are thermalized each- 
other through a well-known mechanism of elastic scattering r e (that defines the conduction 
electron reservoir), at the femto-second time scales (or below), and are also contacted to the 
lattice through the Fermi-Dirac distribution, and inelastic scattering r p h. On the other hand, 
the ferromagnetic order parameter is contacted to the lattice with a well-known relaxation 
time r that is measured in ferromagnetic resonance (FMR) experiments, and is typically of 
the order of the nanosecond (or few hundreds of picoseconds). This description leads to the 
model depicted in Fig 1(b). 

The basic idea developed below lays on the fact that the typical time scales of the dynam- 
ics of the two sub-systems are largely separated. There is a slow variable, the magnetization, 
and fast variables, the degree of freedom related to the spin of the conduction electrons. It 
is then possible to reduce the action of the fast variable to the role of an environment with 
regard to the magnetization, like for spin-bath relaxation. The effect of the coupling to the 
spin-dependent electronic sub-system will then be reduced to specific damping and fluctua- 
tion terms added to the usual stochastic equations for the magnetization. This will be our 
line of reasoning followed in the last section, after describing the two sub-systems. 
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FIG. 1: Thermokinetic picture of irreversible spin-transfer. Ferromagnetic system (with mag- 
netization M), and electric system with spin accumulation density An and electronic density at 
the Fermi level no- The chemical potential /x is defined for each spin channel. The three sub- 
systems are coupled together through the relaxation times t S( { (interband s — d like relaxation) and 
r s f . (intraband spin-flip relaxation). The sub-systems are also coupled to the current generator I, 
and to the heat reservoirs, through the corresponding well known relaxation times tq (Neel-Brown 
waiting time), r e and T p h- elastic and inelastic electronic relaxation times. 

II. SPIN-DEPENDENT TRANSPORT 



In order to explain the high resistance and the high thermoelectric power observed in 
transition metals, Mott introduced the concept of spin-polarized current and sug gest ed that 
s-d interband scattering plays an essential role in the conduction properties [52 1. This 



approach in terms of two conduction bands 0], explained the existence of a spin-polarized 
current in the 3d ferromagnetic materials [53], and was used for the description of anisotropic 

nfi n 

magnetoresistance (AMR) [5_4j,|5Ja], the description of spin-polarizer [56], and thermoelectric 
)ower (5?]]. With the discovery of giant magnetoresistance (GMR) [58^ and related effects 



(like domain wall scattering |6CJ, |6l|, |62j, |63j discussed below ), the development of 



spintronics focused the discussion on _jqpin-fhp scattering occurring between spin-polarized 

conducting channels 0,0,0,000- The two-channel model, which describes the 
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conduction electrons with majority and 
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mmoriw 
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73, 



spins, is applied with great efficiency to 



74j,|75(, including metal/semiconductor 



GMR and spin injection effects 

[3| and metal/superconductor interfaces |7jj. In this context, it is sufficient to describe the 
diffusion process in terms of spin-flip scattering without the need to invoke interband s-d 
scattering. 

It is convenient to generalize the two spin channel approach to any relevant transport 
channels, i.e. to any distinguishable electron populations a and 7 (defined by an internal 
degree of freedom). The local out-of-equilibrium state near the junction is then described 
by a non-vanishing chemical-potential difference between these two populations: A/i a7 = 
/Xq, — /i 7 7^ 0. In other words, assuming that the presence of a junction induces a deviation 
from the local equilibrium, the a and 7 populations can be defined by the a — > 7 relaxation 
mechanism itself, that allows the local equilibrium to be recovered in the bulk material 
(lim z ^± 00 A^(z) = 0) [42]. Such considerations have been presented in some im por tan ; 

Hintronics studies on the basis of microscopic calculations 0, y, y, y, y, y, b y, 
, 0]. The thermokinetic approach [81] allows us to deal with interband relaxation on 
an equal footing with spin-flip relaxation, with the help of the transport coefficients only. 
For this purpose, the two spin-channel model is generalized, with the introduction of the 
corresponding transport coefficients: the conductivities <j a and 07 of each channel define 
the total conductivity a t = a a + a 7 and the conductivity asymmetry f3 = (cr a — a y )/at, 
the relaxation between both channels is described by the parameter L (or equivalently, the 
relevant relaxation times T 7 ^ a ). 



A. The generalized two channel model 

In the framework of the two conducting-channel model, which includes relaxation from 
one channel to the other, it easy to follow step by step the method described in the first 
section. The conservation laws write (assuming a ID space variable z): 



dn a 


dJ a 


dt 


dz 


<9n 7 


dJj 


dt 


dz 



+ ^ 



(9) 



where n a and n 7 are the densities of particles in the channels {0,7}. 
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FIG. 2: Two channel model, including relaxation that couples the two electronic populations. 



The entropy variation writes: 



TI 



'J a 



q ^7 q ^cry^a ^7 ) 



(10) 



the application of the second law o 
coefficients a a > 0, cr 7 > , and L > 



thermodynamics leads to introduce the Onsager 



42 



such that: 



a a dn c 



e dz 
e <9z 



(11) 



where ^ a7 describes the relaxation from the channel a to the other channel 7 in terms of 
velocity of the reaction a — > 7. It is not necessary, in what follows, to distinguish between 
the electric part and the pure chemical part of the electro-chemical potentials (see 82| ). The 
effects of the electric charge distribution are described in Appendix A, with the introduction 
of the screening length I and the relation to the relaxation times. As shown in Appendix A, 



the Onsager coefficient L is inversely proportional to the electronic relaxation times r, 



a<— >7 



L oc 
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f 



(12) 

' cr^7 ' 7— *a J 

where / and g are two functions close to unity, and related to the electric charge distri- 
butions (see Appendix A). Note that due to our definition of \i a and /i 7 , there is no direct 
coupling between the two channels : there is no transport coefficients that couples the two 
first equations in Eq. (fTTj) . This is a consequence of the definition of the electronic popu- 
lations, through the relaxation process itself (the populations are stable if \I> = 0). Indeed, 
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FIG. 3: Junction between to layers / and II. Chemical potential profile over the interval [A, B] 
in the a and 7 channels. The A and B points verify /x a (A) = /Xy(A) and /x a (B) = fx-y(B). The two 
straight lines represent the $ variation in each region (<3?/, &u). It can be directly seen that the 
out-of-equilibrium resistance R ne is determined by the $ discontinuity at the interface. 

the out-of-equilibrium configuration at the interface is quantified by the chemical affinity 
A fx = Li a — /x 7 , i.e. the chemical potential difference of the reaction. 
The total current J t is constant: 



However, it is not possible to measure separately the different conduction channels, since 
any realistic electric contact short-cuts the two channels. What is measured is necessarily 
the usual Ohm's law, J t = —cr t ^, that imposes the reference electric potential $ to be 
introduced, together with the total conductivity o t = a^ + a^ [83]. The potential $ is hence: 

e<3> = — (a a Lt a + a 7 yU 7 ) (14) 

Let us assume that the two channels collapse to a unique conduction channel for a specific 
configuration, the reference, which is a local equilibrium situation: Afi eq = 0. The out-of- 
equilibrium contribution to the resistance, R ne , is calculated through the relation: 




(13) 



J t eR 




— (/i 7 - e$(z))dz 



(15) 
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R ne = ~ I B aa ~^ d ^ dz (16) 



so that 

R ne = - 

J t e J A 2a t dz 

where the measurement points A and B are located far enough from the interface (inside 
the bulk) so that A/a(A) = Afi(B) = (see Fig. EJ). The integral in Eqs. (115)1 is performed 
over the regular part of the function only ($ and <7j are discontinuous) 84]. Eq. ITHl allows 
the out-of-equilibrium resistance at a simple junction between two layers (composed by the 
layers I and II) to be easily calculated. If the junction is set at z = and the conductivities 
are respectively a\ and a\ l (i = {0,7}), we have: 

jTeR ~ = f'AzA^ dz+ [°°J^laA^ dz 

J a 2(7 t dz J 2a t dz 

The equilibrium is recovered in the bulk, so that: 

K»=[- a yT-^-/ jI )— ( 18 ) 

The chemical potential difference Afi(z), which accounts for the pumping force opposed 
to the relaxation a — > 7, is obtained by solving the diffusion equation deduced from Eqs. 
(fTT|) and Q, and assuming a stationary regime for each channels, = ^ = 
HflQ: 




42. 



70 



dz 2 I 2 { ' 

° Z L diff 

where 

l d ? Sf = zL{a- 1 + a- 1 ) (20) 

is the diffusion length related to the a — > 7 relaxation. 

At the interface (z = 0), the continuity of the currents for each channel writes </„(0) = 
Ji'(0), were 

U0) = -^^ + ^J t (21) 

ea t dz a t 

which leads to the general relation: 



ai o*J\ ( oio* a" a^ 1 



-1 



M0)H5-5l h^ + iif I eJ t (22) 



°> a + I \ a t L diff °t L diff 
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Inserting Eq. (|22J1 into Eq. (jl8j) . we obtain the general expression for the out-of- 
equilibrium resistance (per unit area) produced by the a —>■ 7 relaxation mechanism at 
a junction: 



R r - 



2a\ 




where we have used the relation 




<j"a"eL" 



a "7 



a 
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(23) 



r 1 

l diff 



eL 



(24) 



a,(l-/3 2 ) 

It is convenient to describe the conductivity asymmetry by a parameter (3 such that 
o~ a = cr* (1 + /?)/2 and a 1 = a t (l— /3)/2. The out-of-equilibrium contribution to the resistance 
then takes the following form: 



R" 



(Pi -P. 
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(25) 



In the case of the subsystem described by two spin-channel, the relaxation leads to a 
spin- accumulation effect A/x^ at the interface of a two identical ferromagnet with antiparallel 
configuration. The corresponding resistance contribution is: 



l ' art 



PI 



01 



o: t (l-Pl)" Sf y/eLa t (l-ft) 
This expression is the well-known giant magnetoresistance contribution 



(26) 



73. 85 



86]. 
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B. The four channel approximation 



In the previous subsections, two different electronic relaxation mechanisms have been 
invoked separately in order to describe giant magnetoresistance or anisotropic magnetore- 
sistance. It is clear however that the two relaxations would take place in parallel, leading to 
a more complex redistribution of spins within the different channels. In the present subsec- 
tion, we consider a system in which the two mechanisms coexist, leading to a four channel 



model 



44] 
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The generic band structure (energy as a function of wave vector k for a given direction) 
of a 3d ferromagnet is schematized in Fig. 0] The band s is parabolic and the exchange 
splitting is very small. In contrast, the d bands are strongly shifted between up and down 
spin carriers. The hybridized zone is schematized by the dotted lines at the intersection. 
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FIG. 4: Generic band structure for a 3d ferromagnet with s and d bands schematized for an 
arbitrary direction of the wave vector k. The shift between the two d bands for the two spin 
carriers up and down is exemplified. The hybridized zone is schematized with dotted lines at the 
junction between s and d bands. At the Fermi level four different electronic populations can be 
identified. 

The system is composed by the reservoirs of the injected s electrons and the ferromag- 
netic layer composed by the d electrons. At the interface, current injection leads to a 
redistribution of the different electronic populations that are governed by spin polarization 
and charge conservation laws. Let us assume that the current injected is spin polarized in 
the down polarization (!). The conservation laws should be written by taking into account 
the reaction mechanisms between the different populations. At short time scales (electronic 
scattering) the relaxation channels are assumed to be the following four 

(I) e s j_ — > e^i (spin-conserved s-d scattering) 
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(II) e s | — > e s ] (spin-flip scattering for the s population) 

(III) e s j — > e^i (spin-flip s-d scattering) 

(IV) e<jj, — > e<jf (spin-flip scattering for the d population) 

Process (I) is assumed to be the main mechanism responsible for anisotropic magnetore- 
sistance (AMR). Process (II) leads to the well-known spin-accumulation effect and was also 
described in detail in the first subsections. According to the fact that the majority-spin d 
band is full and lies at a sizable energy below the Fermi level, the current J^t is negligible and 
the channel d f is frozen. Processes (III) and (IV) are hence negligible 56]. Consequently, 
we are dealing with a three-channel model {s |, s l,d j}. 

The total current J t is composed by the three currents for each channel : J t = J s f + 
J si + Jdi- m order to write the conservation laws, the relaxation rate ^ s di is introduced 
to account for s — d spin-conserved scattering, and the relaxation rate \l/ s , is introduced in 
order to account for spin-flip scattering. Assuming that all channels are in a steady state 
(this condition will relax in the last section, where the magnetic system is coupled to the 
channels d {) : 



dn t 


dJ t 


dt 


dz 


<9n sT 


dJs^ 


dt 


dz 


dn sl 


dJsi 


dt 


dz 


dridi 


9Jdi ^ 


dt 


dz 







= 



(27) 



sd 



' sd 







where n t , n s i, ra s |, rid[ are respectively the total densities of particles and the density of 
particles in the in the s f, s |, d J, channels. The system is described by the number 
of electrons present in each channel at a given time, that defines the four currents, plus 
the entropy of the system. The conjugate (intensive) variables are the chemical potentials 
{/i s -f, fj, s ±, fidh Hdi}- As described in Appendix B, the application of the first and second laws 
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of thermodynamics allows us to deduce the Onsager relations of the system : 



Jsl 
^sd 



si 



e dz 

_avr fyvr 

e dz 
°di d^i 
e dz 

L s d (Vsi - f*dl) 



(28) 



L s (/^sT ~~ Mai) 

where the conductivity of each channel , <7 S |, cr^, cr^} has been introduced. The 
first four equations are nothing but Ohm's law applied to each channel, and the two last 
equations introduce new Onsager transport coefficients (see Appendix B), L sd i and L s , that 
respectively describe the s — d relaxation (I) for minority spins under the action of the 
chemical potential difference A/ij_ = /x s j./2 — ^i and the spin-flip relaxation (II) under 
spin pumping A^ s = fi s ^ — fi s i/2. According to Appendix A, the Onsager coefficients are 
proportional to the corresponding relaxation times. 

For convenience, we define the usual charge current J 0s = J s \ + J si-, the minority-spin 
current J | = Jsl + Jdii an d the two polarized currents 5J± = J s ± — J^i and 5J S = J s ^ — J s ^. 
We introduce the o~ s and o~-\ conductivities {cr s = a s ^ + a s i and <7j = a s i + o"dj}. The 
conductivity imbalance /3; and f3 s between respectively the s j and d j channels and the s | 
and s I channels are: 



Eqs. (|2Tjl becomes : 



Pi 

Ps 



dj_ 

dz 
dJoi 
dz 
dSJi 
dz 
dJ 0s 

dz 
d5J s 

dz 



Q~sl - Q~dl 

0~ R 



dJ dl . dJ x 



(29) 



+ 



dz dz 

i>s 

-2tp sd - ip s 

-1psd 

ipsd - 2^ s 



(30) 
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and, denning the quasi-chemical potentials fi s = /i s j + /i s |/2 and /x^ = n s j2 + /i^j., Eqs. f)28|l 
becomes : 



Jqs 
^sd 



2e\dz Pl dz 
°l (g 9 H , ^A/i^ 

2e V 0* <9^ 



<9z 



2e V <9z 



?± (r ^tl _|_ 

2e I s 9z <9z 



(31) 



The equations of conservation [Eqs. ([3(jp] and the above Onsager equations lead to the two 
coupled diffusion equations : 



where 



dz 2 



1 



\ 2 



(32) 



A solution of Eqs. 



A, EE 



hf = 



W(l-lf) 

f ^(l + /3|) 
2eL s 

4eL s 
'cr a (l-ft) 



2eL 



6(1 



(33) 



is 




A/ii + A/i 2 



A! I I IT - TT A// 



/ 2 

<sd 



(34) 



19 



with 

/ z z_ 

= aie A + + a 2 e Aj * 




(35) 

&ie A - + b 2 e A - 



where 

_ , , __ , : ^ I 1 -U 1 A 

2 ^sd ^ l s/ 



Ai 2 =k? + tf)|l± 



I — I — -2 



^ te? +',-/)• 

The constants a%, a 2 , b\, b 2 are defined by the boundary conditions. It can then be seen 
that the usual spin accumulation corresponding to Afi s also depends on the spin-conserved 
s — d electronic diffusion which is known to be efficient and, conversely, that spin- 
conserved diffusion is able to lead to a spin accumulation, or d spin- accumulation effects. 
Accordingly, we expect to measure some typical effects related to spin-accumulation in single 
magnetic layers, or if (3 S = : this point will be illustrated in the new expression of the 
magnetoresistance (Eq. (J3*§|) below), and in Section IV through the effect of current induced 
magnetization switching (CIMS). s — d relaxation adds a new contribution to the resistance, 
which plays the role of an interface resistance arising from the diffusive treatment of the 
band mismatch 

The resistance produced by the usual spin-accumulation contribution, plus the contribu- 
tion of s — d relaxation, are defined (see Eq. (fIH)l) by 

Rsa = ~^ [ |- {JM ~ dz (36) 
eJ t J B dz 

where $>(z) is the total electric field and fa is one of the chemical potentials. Providing 
that the total current is J t = J s ^ + J si + Jdi, or 

Jt = ~- — fidi + — Vsi + — /V 37 

e dz \a t a t a t ) 

The total electric field can also be written (from Eqs. (|28|)) as 

a t e \ a t dz a t dz a t dz ) 
where a t = a s ^ + cr s j_ + o~di- The resistance is given by : 

R, a= -\ [ B (^ a ^i + ^ a A^.) dz (39 ) 
eJ t J A V °t dz o t dz ) 

This three-channel model brings to light the interplay between band mismatch effects and 
spin accumulation, in a diffusive approach. It is interesting to note that the local neutrality 
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charge condition which is often used (see for instance Eq. (4) in |87[) was not included, as 
described in Appendix A. On the contrary, we have imposed the conservation of the current 
at any point of the conductor. Indeed, electron transfer from a channel to another where 
the electron mobility is different, induces a local variation of the total current. 
The resolution of the coupled diffusion equations is discussed elsewhere j^J. 

C. Domain wall scattering 

In the description performed until now, the spin quantification axis that defines up and 
down spin states was fixed through the whole structure (i.e. through the layers and the 
interfaces). Providing that the spin quantification axis follows the direction of the magneti- 
zation, it could be non-uniform throughout a ferromagnetic layer, or crossing an interface. 
This is especially the case in the presence of a magnetic domain wall. In a thin enough 
magnetic domain wall the spin would not follow adiabatically the quantification axis, lead- 
ing to spin-dependent domain wall scattering (DWS) jfifl Ifill l^. . This effect has been 
investigated intensively in the last decades in various structures 69]. The underlaying idea 
is however rather simple, and can be formulated easily with a generalization of the two-spin 
channel approach. For the sake of simplicity, this generalization will be performed only for 
the two electronic populations {a, 7}. 

As performed in reference |42j (and appendix B), we start with the conservation of the 
particles for the two channels, in a discreet model. The system is described by a layer in 
contact with a left layer and a right layer The spin-flip scattering introduced in 

the previous sections is described be the reaction rate A probability (1 — Ae(k)) of spin- 
flip alignment along the quantification axis is introduced. In the case of ballistic alignment 
(1 — Ae(k)) = cos 2 (A6(k)/2) where A6(k) is the angle between the magnetization of two 
adjacent layers and £&. The conservation of the particles is now describes by: 

^ = (1 - Ae(k))I k a -^ k - I k ^ k+1 + Ae(k)I^ k - V k 

dN • ^ 

^ = (1 - Ae(k))I k -^ k - I k ^ h+1 + Ae(k)I k -^ k + * k 

With the notation introduced in the previous sections, the entropy variation can be 
written in the following way (Appendix B): 
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T 



dS_ 
~dt 



r>R l ->l pf2->iT 



n 1 

5 (A/" 1 - A/ + 2(1 - Ae(*)) V) tf/*" 1 " 



k=2 

n 



(41) 



fc=2 



fc=l 



where we have introduced Jo 
and R - 1 



Jq, + Zy, 51 — I a — Ij, fiQ 



The terms P, 



. i/j. — j, ; i. - . /Mi — /^q "I - /x-y, and A/i — /Xq, zi^,. 
stand for heat and chemical transfer from the reservoirs to the 



system E. 

After performing the continuum limit, the internal entropy production X (or irreversibil- 
ity) reads: 



T.l 



1 dfio 1 

2 97 Jo+ 2 



dz 



+ 2eA/i 5J + AfiV 



(42) 



The first term is the Joule effect, the second is the dissipation terms related to the spin- 
accumulation process that occurs at the interface, or for magnetic domain wall, and the 
third term is the dissipation due to spin-flip (or s-d) electronic relaxation. The expression 
of the entropy production Eq. (}4*2*j) allows the Onsager relations generalizing Eq. ([lip or 
Eq. H2SI to be deduced: 



Jo — — — 



5J 



2e 

<J0 



dz 



2e V oz oz 



(43) 



where e = HfflA^o^) and, as already introduced: = a a + <J y and (3 = (a a — ap)/a t - 
The diffusion equation for A/i, obtained in the stationary regime, is modified accordingly: 



<9 2 A/i 



1 



' nw / 



A/i 



IdA/i 



where the length l di ff as been defined in the first section Eq. (|2Uj) 



(44) 



i. 



^(l-/? 2 ) 



2eL 



(45) 
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the domain wall diffusion length Idw is defined as: 



while the length k is given by: 

2(3 2 

The magnetoresistance is modified with respect to Eq. ([19|h due to the new term <9 2 A/x(z) 
in the diffusion equation. It is worth pointing out that a spin accumulation Afi(z) ^ est 
should be expected in case of spin polarized current {(3 7^ 0) even without the usual spin-flip 
contribution, i.e. in the ballistic limit. 



III. FERROMAGNETIC BROWNIAN MOTION AND MAGNETIZATION 
SWITCHING 

A. Thermokinetic derivation of the Fokker-Planck equation 

The description performed in the previous sections is related to the transport properties 
of charge carriers in case of spin polarized current. In spintronics experiments, the electric 
current is spin-polarized through a ferromagnetic layer, but it is not necessary to describe 
the ferromagnetic order parameter as such. This is of course no longer the case for current 
induced magnetization switching experiments, where the magnetization is the measured 
variable. 

The magnetization is a fascinating degree of freedom, that has to be described in length 
in terms of rotational Brownian motion. The description of the dynamics of ferromagnetic 
particles coupled to a heat bath is a very active field of investigation 



des of ferromagnetic 

00000, 



and the resulting predictions are rat 



rer well known and validated experimentally at large 



000 and short time scales. The magnetization relaxation 

described here is limited to the so-called Neel relaxation that involves only the magnetic 
moment, in contrast to the Debye inertial relaxation occurring in ferrofluids (in which the 
ferromagnetic particles rotates in a viscous environment, leading to surprising inertial effects 
like negative viscosity 0]). 

The aim of this subsection is first to show that the rotational Fokker-Planck equation 
governing the dynamics of the magnetization M of one monodomain particle coupled to the 
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heat bath can also be obtained applying step by step the approach used in the previous 
sections. The resulting Fokker-Planck equation with the corresponding Onsager transport 
coefficients, and the hypothesis performed, can then be compared term by term to the 
previous study of spin-dependent charge transport. 

1. Geometrical representation of the statistical ensemble 

Let E be a statistical ensemble of N identical monodomain particles of volume v, having 
the same energy per unit volume V ma9 (8, 0), magnetization M and thermostat temperature 
T. The vector M is defined by the angles 9 and <fi. The ensemble £ can be represented by 
a distribution of representative points over the unit sphere (fig. Ej) with a density n(9, </>). 




FIG. 5: a) The figure from the left illustrates the flow of representative points over the unit sphere: 
J and J^. b) The figure from the right illustrates a particular case of distribution of points on 
the sphere: the points are concentrated at two attractors, one with more particles than the other 
(asymmetric double well potential). 

We divide the ensemble of representative points £ in sub-ensembles T, e ^ such that the 
magnetization is confined within the solid angle dVg^ = sinOdOdcj). (i.e. the representative 
points lie between two consecutive parallels and meridians over the sphere). 

As the particles undergo changes of magnetization orientation, the representative points 
move on the sphere, and there is a net surface flux of representative points J ma #; the repre- 
sentative points move from one sub-ensemble to another sub-ensemble ^e+Ae,4>+A4>- The 
probability of finding a particle with the magnetization orientation within the solid angle 
dVerf at a given time t is dP(6, 0, t) = "^'^ dVe^. 
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2. Conservation laws 



The sub-ensembles of representative points ^e,<f> are described by the following extensive 
parameters: the entropy dS = s(9, 0, t)dVg^, the number of points dN = n(9, </>, t)dVg^ and 
the energy dE = e(9,(p,t)dVg t( j ) , where s and e are the entropy and energy densities. The 
flow (of points, energy, and entropy) is described by the flux J (J n , J e and J s ): 



J=J% + J% (48) 

and accounts for the flow of the corresponding magnetic moments relaxing or precessing 
along the coordinates 9, 0, where u r , ug, are the unit vectors in the spherical coordinate 
system. 

The conservation laws of the number, energy and entropy of the particles contained in 
the sub-ensemble S e $ write : 



f dn 
~dt 



div J„ 
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— = — div J e 
dt 
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— = - div J s + 1 
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(J e n sin9) ■ 

(J e e sm9) 

{J e s sin9) 



1 dJt 



sin9 
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dJ* 


sin9 


dcj) 


1 


dJ* 


sin9 dd> 



(49) 



+ 1 



where in contrast to the energy and number of particles, the entropy s is not a con- 
servative quantity, and an internal entropy production term I (or irreversibility) is added 
to the entropy flux J s (third equation in Eq. (J49|) ) . 

The expression of the first law of thermodynamics, allows the energy variation to be 
expressed as a function of the partial derivatives that define the chemical potentials and the 
temperature \x = and T = |j . The intensive variables fx and T are also functions of 

(9, <f), t) except if imposed by a reservoir. 

ds 1 de pL dn , . 

di = Tdt~f~dt ^ ' 

where ft = d n e contains all contributions to the energy (see below)). 

The expression of the internal entropy variation can be obtained using the conservation 
laws: 



ds 1,-7 A* ,. ? 
— = — — div J e + — div J n 
dt T T 



(51) 
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or, in an other form: 



^ = -div \ 7j,Je-j, Jn) + Je ■ grad (JjA - J n ■ grad (j^J (52) 

Comparing this last equation with the third eq. from (j49j) . we can deduce the form of 
the entropy production X: 



°S rp °e rj-, 'Jn 



IN f - ,/y3\ < 53 > 



X = J e • I - J - J n • ^— 

The entropy production X is a sum of products between the fluxes Jk and the correspond- 
ing conjugate forces Ff. 0|. 

We assume in the following that the temperature T(9, 0) = T is fixed by a unique 
thermostat: the first term in the right hand side of the eq. vanishes. 

A sufficient condition to impose the second law of thermodynamics X > is then to 
build a quadratic form. This leads us to define the matrix £ of Onsager transport co- 
efficients Lij(6,(p) (that are state functions of dimension [energy]~ 1 [time]~~ 1 ) such that 
Ji = £j- (Lijdjjl), where the symmetrized £ matrix is positive. The Onsager reciprocity 
relations impose furthermore that = ±Lji, where the sign (-) is present if is a func- 
tion of the magnetic field (there is no angular velocity here) [35]. 

The following relations are deduced: 

{ e _ dfl 1 dp, 

Jn -- Lee de- Le ^d^ 
, on i op (54) 

where Lq^ = —L^e- The first of equations (j4Uj) re- writes: 



dt 



-div J n = +div (L gradji^j (55) 
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(56) 
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where Lee > and > 0. This is the general expression of the density variation 
dtn(8, <p, t) of particles number from the sub-ensemble E^. Note that there is no relaxation 
terms in the conservation law (J55|) of the magnetization: the flow of representative 
points is conserved on the unit sphere. This assumption will be removed in the case of spin 
injection performed with electric currents (see next section below). 

In the same manner as for Eq. (j!2|) . the Onsager coefficients are related to the relevant 
relaxati on ti mes. In the case of ferromagnetic insulators, the relaxation channels are well 



defined 



10l| . and the coefficients L^ are directly related to the relaxation times T\ and T 2 



measured in ferromagnetic resonance experiments. 



B. The rotational Fokker-Planck equation 

In thermokinetics, the intensive parameter which controls the number of particles of a 
sub-ensemble is the chemical potential ft. The relevant energy terms are contained in the 
deterministic potential /x, and the stochastic term is defined by thermal fluctuations due to 
the coupling to a relevant heat bath. Anticipating the last section, it is worth pointing out 
that the relevant heat bath is defined by the degrees of freedom of the environment which 
are that of the lattice or that of the electronic system (as discussed below). The fluctuations 
are taken into account through a temperature dependent chemical potential that takes the 
following form (derived in the general case by P. Mazur in Ref. jssj) : 

fl = k B Thx^+vV ma9 {d,(f>)+^ (57) 

The first term in the right hand side of Eq. (JSTj) is responsible for thermal agitation at 
temperature T, the second term vV ma9 represents the magnetic energy of one particle that 
defines the local magnetic field H e ff = —d^{V ma9 ) and the third term is a constant which 
is related to the chemical nature of the particles. 

The local equilibrium condition dip, = defines stationary flux (due to both drift and dif- 
fusion) that are mutually compensated along the coordinate i. This point is well illustrated 

n 

in the work of Guggenheim while introducing the electro-chemical potential [33| in order to 
generalize the description of an electric fluid to ionic solutions. 

Inserting the expression of ft, and using the reciprocal relation Lg^ = —L^, the equations 
for fluxes and the variation of particles take the form: 
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/ Qymag g , Qymag \ fk B T8n g' k B T 8n\ 
\ 89 sinO d(j) J \ v 89 sinQ v 8<f) J 

/ /d ymag ft y Qymag \ / fa? fa ft 8u\ ^ ' 

\ 9 89 + sine k B T 8<p ) " ~~ K 9 V¥ + sineM ) 
where the following standard notations have been introduced: 

h' = ^>0; j = -h£L = h*!L. k' = L^>0 (59) 

n n n n 

Assuming that g' is constant, Equation Eq. (|56j) rewrites: 
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(60) 



The expression Eq. (I60|) represents the rotational Fokker-Planck equation obtained by 



thermokinetic means. Its expression is identical to that obtained by Brown [89j through 
stochastic calculations. 

Furthermore, the Onsager matrix also follows the symmetry of the system, and is invariant 
by rotation around the anisotropy axis, so that = Le,e- The following relations are 

obtained : 

Lap = Lj» 
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k B T 



(61) 



which permits us to write the equations above in a compact vectorial form: 
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(62) 



where the gradient V is in spherical coordinates and u r is the spherical radial unit vector. 
It is to be noticed that the second equation has drift terms which contain VV ma9 , and 

diffusion terms which contain Vn. The terms k' = h'^- and q'^- are the rotational 
diffusion coefficients and the terms g'n and bin represent drift coefficients. 
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C. Landau-Lifshitz-Gilbert equation with diffusion 



Furthermore, using the first equation from (|62p. one can deduce the Landau-Lifshitz- 



Gilbert equation with diffusion. As J 



n 



du, 
dt 



, we arrive at the equation: 



du. 



dt \ v n 



u r x 



r >< + 



v n 



(63) 



where the first term in the right hand side is the precession term, and the second term in 
the right hand side describes the longitudinal relaxation. Multiplied by the amplitude of the 
magnetization M s , becomes the Landau-Lifshitz-Gilbert equation with diffusive terms: 



g'M s M x H eff - -=h ) + ti 



M x I H eff - 



ksT Vn 

vM a n 



x M 



(64) 



dt " * \ 11 vM s n 
Experimentally, the first contribution can be observed through ferromagnetic resonance 
measurements (FMR) at typical frequencies of tens of GHz (around 100 psec time scales). 
The thermalization time (proportional to 1/g' ; see next section) is given by the width of 
the resonance peaks. Both frequency and time resolved noise experiments have been also 

99,0. In 



96, 



97 



98, 



performed in order to measure precession and thermal spin-waves 
the above measurements, the data are averaged over many shots (or trajectories) near an 
equilibrium position (linear response regime). In contrast, the measurements at large time 
scales, typically beyond few nanosecond, access the magnetization reversal for which the 
precession terms can be neglected. The one shot measurements give a direct access to the 
stochastic nature of the signal 21, 3, 25, a snapshot is a statistical event, namely 

the magnetization reversal from one metastable state to the other, that is governed by the 
random fluctuations, described by a "Langevin force" that is not present in the averaged 
LLG equation. 

Equation Eq. (jo3|l can be put into the Gilbert form by performing the cross product xM 
at the left and right hand side of the equation. We obtain the well-known Gilbert equation, 
that defines the Gilbert damping parameter rj: 



dM ^ „-> / -* dM 
TMx [ H eff -r)- 



dt 



dt 



(65) 
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where T is the gyromagnetic factor. The constant h', g' and k' are related to the Gilbert 
damping coefficient rj, T and the magnetization at saturation M s through the following 
relations: 



a (1 + ( V TM s y) M & 

or 



> 0; h' = r]TM s g f ; (66) 



^fk\jm-^° (67) 

D. Activation regime and Neel-Brown law 

1. Neglecting precession 

In the slow relaxation measurements (the so called magnetic after-effect), relaxation is 
governed by activation over a potential barrier. At longtime scales (beyond tens of nanosec- 
ond to hour), the precessional terms can be neglected. The expression for the surface current 
fluxes Eq. (J58|) becomes 

Qymag i ,k B T 8n 

DO V 09 (a.o\ 

Qymag f k B T On 1 J 

and the corresponding Fokker-Planck equation becomes: 
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(69) 



These expressions will be used in the following paragraph for deriving the Neel-Brown 
relaxation time. 



2. The double-well potential and the relaxation times 

The double well potential (see fig. EJ) is the first approximation of the ferromagnetic 
particle energy V mag (9, 0) beyond the harmonic potential, but it is also a realistic magnetic 
potential in the case of an uniform magnetization with uniaxial anisotropy [2^ : 

ymag = R ^2 g _ cos (fl _ ^ ( 70 ) 
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where 9 is the direction of the magnetization, ip the direction of the applied field H ( 
and K is the anisotropy constant in energy per unit volume. 
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FIG. 6: Double well potential (continuous line) with stochastic fluctuations (dashed area), and 
the definition of the barrier heights. 

In order to apply Brown's method, we consider a potential which has minima at 9\ and 



9 2 = ~K—9x, and a maximum at 9 m . Following Kramers transition theory, Brown 89J assumes 
that most of the representative points on the unit sphere are concentrated at the energy 
minima of V mag {9) where they are in thermal equilibrium so that locally n takes the form 
of the Maxwell-Boltzmann distribution. Thus only a minute fraction of the representative 
points is outside the energy minima allowing a small diffusion current between them so 
manifesting the non-equilibrium conditions. 

{ n(9 1 )e-^ [V(e) ' V{ei) \ for 9 E (9, - e, X + e) 
n(0,<f>) = { (71) 
{ n(9 2 )e'^ [V(e) ~ Vie2) \ for 9 e (9 2 - e, 9 2 + e) 

The number of particles N\ «-> N 2 from the first well, respectively the second is: 



N x = 2im(9 1 )e^ TV(ei) I u where h = e~^ V{6) sin9d9 



n+e 



n-e 

N 2 = 2im(9 2 )e^ TV{d2) I 2 , where I 2 = e~^ V{e) sin9d9 

Je 2 -e 



(72) 



Assuming that the flow between the two minima Q\ and 9 2 is quasistationary, approxi- 
mated by a divergenceless current 8^1, the total current of particles over the potential barrier 
can be written: 

I = 2nsin9J° (73) 
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Rewriting the first equation fromlHEl one obtains: 



dn v dV ma9 Iv 

d6 + kf dd U ~ ~ 2ixh'k B Tsin6 ^ ' 



which defines the activation regime. Introducing I m as 



Eq. (JZH) yields 



which has the form 



with 



I ^L.,^r de (75) 



I = -S l = & 1 = J»f(*>V) (76) 

Vl m \I 2 h J 



No iVi 

Nl = -N 2 = ^ - (77) 

7~2 7"1 



_ IJ m V 

Tl ~ h'k B T 

hlrnV 

T 2 



(78) 



h'k B T 

Because of the rapid decrease of the exponential factor with distance from the minima of 
V mag , we may in 1%, I 2 , I m replace V ma9 (8) by its Taylors's series about 61,62, respectively 
6 m truncated at the 9 2 term, and replace the upper limit of the integrals by 00. With these 
approximations , we find 

v(v{e m )-v(0 2 )) 

(79) 

t 2 = t 02 e k B T 
where the waiting times are given by the expressions: 

Equation Eq. (J75j) is a formula for a symmetric bistable potential which has minima in 
6\, 6 2 = 71 — 61 and a maximum in 6 m = tt/2. 

For (p 7^ 0, the potential V mag (6, 0) has an as ym metric bistable form, and all arguments 
leading to Eq. (|79j) also apply for an arbitrary 99 j9jj . The general equations are very similar 
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to Eq. (|79|). the only difference in the analytic expression is that instead of the symmetric 
angles Q\, 62, we have the asymmetric angles 6 a, respectively 6b- 

It has to be emphasized that T\ and r 2 are the relaxation times related to the first and 
second potential barrier (starting from the first or the second minima), and that Eq. (J79|) 
constitutes the Neel-Brown law for the particular case of the symmetric bistable well. Inter- 
estingly, in all cases [9]] , the Gilbert damping is reduced to the prefactor of the exponential, 
and consequently plays only a negligible role in the activation process. 

In many cases, the potential is highly asymmetric, and the Neel-Brown law is written in 
the following asymptotic form: 

^.expW^TO (81) 



with three phenomenological parameters a ~ 3/2, AVo and H® w 

0Q- The laws dznj) 

and ()81j) are well established experimentally in usual magnetic sub-micro structures. More 
surprisingly, they have been also observed in CIMS experiments under spin-injection with 
high effective temperature T eff (2000K to 20000K) instead of 300K to 340 K 
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22 



23, 



24, 25] (see next section). 

Fig. describes how to measure the ferromagnetic potential landscape with the help of 
slow magnetic relaxation measurements (the so called magnetic after-effect measurements), 
performed on a single magnetic domain nanostructure [zj. In the example shown, the angle 
and amplitude of the magnetic field is set in order to obtained the two-level fluctuation effect 
(measurements reported in Wernsdorfer et al. 0))- The principle of the measurements is 
sketched in Fig. [7|for a sample with uniaxial anisotropy : the hysteresis loop describes the 
succession of equilibrium magnetization configurations (or quasi-static states) as a function 
of the magnetic field (the field is normalized to the anisotropy) for different angles ip of the 
applied field. The angle 6 describes the direction of the magnetization. For a given angle 
of the field if, the hysteresis is composed of the reversible configurations and two symmet- 
ric irreversible jumps over the potential barrier. The jump occurs from one equilibrium 
(metastable) state defined by the angle 61 to the other, defined by the angle 6 2 (Fig. Efb)). 
At zero Kelvin (without fluctuation), this angle is related to the applied field through a 
relation that depends on the switching mode (the switching mode concerns the magnetic 
configurations occur ring during the jump, i.e. the modes with a typical lifetime of fraction 



of nanoseconds 
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IO2J). In contrast, the amplitude of the jump is easily measured (it 
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is given by the difference M s (cos92 — cos9i)), and gives direct access, from the hysteresis 
loop, to the quasi-static configurations (i.e. to the position of the potential wells). Chang- 
ing the field (amplitude and angle) necessarily changes the equilibrium position (defined 
by W(H) = 0) and the angles 9\ and 62- Note that the case of uniaxial symmetry with 
if = is the unique pathological case, where the change in the amplitude of the field does 
not change the equilibrium magnetic configurations. However, since this pathological case 
is unfortunately that used in most calculations for the sake of simplicity (this was also the 
case here for the derivation of the Neel-Brown's law), the fact the initial and final states are 
necessarily modified is not mentioned in many reports of ferromagnetic relaxation experi- 
ments. Inversely, a change observed in the two equilibrium configurations implies a change 
of the potential landscape, i.e. a change of the effective field. 

The possibility of distinguishing between the effect of an effective field and the action of 
the environment will be of fundamental importance with regard to the study of magneti- 
zation reversal due to spin injection (that are performed on single domain magnetic nanos- 
tructures), because the main problem is to identify the typical relaxation times that govern 
the mechanisms responsible for the magnetization reversal (effective fields vs. stochastic 
fluctuations). It is worth pointing out that it is not relevant (and sometimes misleading) to 
measure relaxation effects in order to access the valley of the potential (i.e. the reversible 
magnetic configurations induced by the action of an effective field) because the position of 
the valley is given by the quasi-static measurements (the hysteresis loop). Furthermore, the 
stochastic fluctuations do not perturb the quasi-static configurations (e.g. reversible part 
of the hysteresis), except at the two critical points where the potential barrier is of the 
same order as the thermal energy. In other terms, only the position of the irreversible jump 
is modified by the activation (or by the observation time window, or equivalently by the 
velocity of the sweeping field [103]), but not the rest of the hysteresis loop. Of course, this 
is no longer the case for larger samples, beyond nanoscopic dimensions, because the sample 
is a distribution of single domain sub-systems with a distribution of anisotropics, defects, 
etc. In contrast, the barrier height - or the amplitude of the thermal fluctuations - can be 
measured only by the activation process, through statistical measurements. If the barrier 
height is above the energy of the lattice kT, the two level fluctuations occur within a typical 
time window that can be tuned with the amplitude of the magnetic field or the temperature. 
In order to access to the relaxation times, statistics should be performed over a significant 
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number of shots. If the exponential relaxation is verified, the relaxation times can then be 
extracted. The ratio of the relaxation times t(#i)/t(#2) gives the asymmetry of the double 
well, and each relaxation time gives the corresponding barrier height. The Neel-Brown law 
is tested by varying the temperature and the magnetic field (Fig. EJ). The whole potential 
can be rebuilt form these measurements. 
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FIG. 7: Slow ferromagnetic relaxation measurements in uniaxial anisotropic single magnetic do- 
main, (a) hysteresis loop at different angles (p of the applied magnetic field. The saturation (i.e. 
initial states without excitation in the relaxation protocol: H y = 0) corresponds to the magne- 
tization Of 11 * or 0| • (b) double well potential after the application of the excitation Aip: the 
whole potential landscape is modified (the barrier height and the two equilibrium configurations), 
(c) two level fluctuations measured between the two configurations (the magnetization jumps from 
potential well at 9\ to the potential well at 62- Right: the Neel-Brown law is verify after perform- 
ing statistics over a significant number of jumps, in order to access to the mean relaxation times 
t(0i) and r(6 2 ). Reprint with the permission from Ref. [9J| Phys. Rev. Lett. 78, 1791 (1997), 
Copyright@American Physical Society 
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IV. SPIN-INJECTION INDUCED FERROMAGNETIC BROWNIAN MOTION 



Let us now consider a system composed of a ferromagnet in which spin-injection (with 
high current density) is performed. In such spin-transfer experiments, the magnetization 
configurations are measured while injecting the current. Usually the magnetic configurations 
are measured with GMR or AMR properties fo r con venience, but micro-magnetometry or 



magneto-optics measurements are also possible [1041 ]. The magnetization is a macroscopic 
(thought nanoscopic) variable described in the previous section (Sec. Ill), and the spin- 
injection is described by the two -or four- channel model presented in section II. We start 
with the assumption that the magnetic system is an open system, composed on one hand 
of the spin-accumulated charge carriers, and on the other hand by the ferromagnetic order 
parameter. It is worth pointing out that in this picture, we are not dealing with a the- 
ory of itinerant ferromagnets, in which spin injection is performed. Such a theory would 
allow the dynamics of the macroscopic variable M to be derived from the spin-dependent 
electronic populations 5n (e.g. a theory in which the LLG equation would be derived from 
the Hamiltonian of the electronic system jioj ]). Such a theory would be very difficult 
because statistical projections should necessarily be performed in order to reduce the micro- 
scop i c degrees of fr eedom to fluctuations and damping for the macroscopic order parameter 
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108, 



109]. The simple phenomenological model proposed here can be viewed as a 
first step in this direction, by defining the different sub-systems, that are coupled through 
relaxation processes (see first Figure in Sec. I). 



A. Current-dependent effective field? 

Before dealing with the dissipative coupling invoked above, a first trivial mechanism that 
couples the two sub-systems should be discussed. This is a deterministic coupling through 
a common electro-magneto-chemical potential /i, from which the associated effective fields 
are derived d x fi , d x Afj, (electric) da[i (magnetic), A/i (spin pumping). This concept follows 
that introduced by Guggenheim j33j when he defined the electro-chemical potential in order 
to describe electrochemical or electrophoretic processes. For instance, the local electric field 
d x Afi is spin- dependent, and might have an action on the magnetic order parameter. On 
the other hand, the magnetic field d^V mag is able to have an action on the spin polarized 
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current. 

These effects are defined by the deterministic part /i of the chemical potential \x and 
can hence be reduced to the action of the effective field H e ff(I) = — V/x (drift part of the 
Fokker-Planck equation, to be opposed to the diffusion part), according to Eq. (J60j) . In other 
terms, within this hypothesis, CIMS experiments reproduce the slow magnetic relaxation (or 
magnetic after-effect) experiments described in the last paragraph of the previous section, 
where the magnetic field excitation is replaced by the injection of current. 

In the context of thermokinetics, the two mechanisms are identical: the current or the 
field drives the magnetization from a stable initial position (# im *) to the metastable state 
Q\ or 82 by modifying the position of the potential landscape. If a relaxation is observed, 
this is due to the usual activation process only, i.e. the jump over the bar rier i nduced by 



the thermal fluctuation (helped by the Joule heating due to the current) HCt 111 lL Ill2 |. 
Accordingly, as underlined at the end of the previous section, this deterministic mechanism 
due to current injection should first be observed on the stable magnetic configurations, by 
measuring the modification of the potential wells (e.g. measuring the modification of the 
hysteresis loop of a single domain nanolayer) under current injection. 

However, it is worth pointing out that the modification of the potential landscape is 



not observed: the basic prediction o 
been verified until now 111 




the determinist action of the spin-injection has not 
Fig. IS] shows two shots with two level 



23, 



fluctuations in a time interval of milliseconds (a) (|2jj) and microseconds (b) (|23j), obtained 
by different groups with GMR measurements on trilayer nanopillars, and single irreversible 
jumps measured in an hysteresis loop of a Ni nanowires (c) (|25j): each point in the zoom 
(right part of Fig. El( c )) is a 6 microseconds pulse. The hysteresis is performed with current 
varying from 2.4 10 7 A/cm 2 to 1.5 10 7 A/cm 2 (the position of the jump without current 
injected is shown by the arrow). In all cases, the metastable states (described by the angles 
8± and 82 defined in the previous section) coincides with 8 1 ™ 1 and 8 l 2 mt measured before the 
application of the current. All happens as if the potential landscape were not modified! 

In contrast, the action of the current is huge (some fraction of an eV to few eV, i.e. from 
25 % to more than 100 % of the anisotropy energy of the ferromagnet), and is observed only 
for the irreversible jump. By acting on the jump and not on the potential landscape, the 
effect of the current mimics the action of a temperature. 

The arguments developped above are based on an analysis performed with neglecting 
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FIG. 8: Ferromagnetic relaxation measurements under spin-injection, (a) Two level fluctuation 
measured with GMR in a trilayer system (Permaloy/Cu/Permaloy, from 21]). (b) Two level 
fluctuations measured with GMR in a trilayer structure (Co/Cu/Co, from 23]). (c) Hysteresis 
loop measured with AMR (left) and zoom around the irreversible jump measured under spin- 
injection (6 microseconds pulse per point) (details reported in 25j|). The hysteresis measured 
with AMR shows the succession of the equilibrium configurations. In all cases, the initial and final 
states (the equilibrium states) are that measured without current injection (9 = 6 imt ). Equilibrium 
configurations are not modified by the spin-injection. Reprint with the permission from Ref. 21] 
S. Uhrazdhin, et al. Phys. Rev. Lett. 91 146803 (2003) and A. Fabian, et al. Phys. Rev. Lett. 91 
257209 (2003), Copyright® American Physical Society. 

the precession terms. If we assume an hypothetical precession maintained in a stationary 
regime, the magnetization is then driven by the precession, and the trace in Fig. 6 should be 
interpreted as trajectories in the phase space. Within this context, the precession induces 
intermittency, i.e. a highly specific chaotic behavior in which the time spent out of the 
attractors is negligible. Such an approach is discussed in some theoretical works about CIMS. 
This interpretation is however evacuated here because the measured traces (Fig. 6) mimic 
exactly the full stochastic process composed of the two-level fluctuation (that follows the 
Neel -Brown activation law), superimposed to the same white noise for both quasi-static states 
1 . The observation of identical noise in both states is in contradiction with published 
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1151] . Furthermore, 



simulations performed with the hypothesis of spin transfer torque 
after spending years in measuring activation processes in ferromagnetic nanostructures, it 
is difficult to accept that a deterministic behavior replaces the activation and mimics so 
perfectly the full stochastic behavior. 

Accordingly, we focus our attention on the dissipative, or irreversible, spin-transfer effect, 
for which the action of the current is expressed through the diffusive terms of the Fokker- 
Planck equation . 



B. Negative damping 




FIG. 9: Under spin-injection, the ferromagnetic layer is an open system coupled to the spin- 
accumulation reservoir. The variation of the density of magnetic moments (oriented at a given 
direction in the unit sphere) is given by the divergence of the magnetic flux added the correction 
due to the relaxation rate \& of spins relaxing from one system to the other at the interfaces. 

In order to describe the ferromagnetic layer coupled to the spin-dependent electric sub- 
system, we rewrite the conservation law of the ferromagnetic systems (Eq. IB3j) with adding 
the relaxation term \l/ due to spin-transfer, while substracting it from the spin-accumulation's 
conservation law: 
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g n f erro 

dt = - dlv(J n ) + * 

() " -div(5J)-^f 



i /( r (^^) 



dt 

According to the first law of thermodynamics, the contribution due to the electronic 
relaxation \I/ should be taken into account in the entropy production of the ferromagnet. 
Inserting Eq. (jH^j) into Eq. (JoTj) . the entropy production is: 



X = J e ■ grad (I) - J n ■ grad (~ \ + 



(83) 



where A/x = /i e ; ec — /2 ( e.g. A/i = // s j — /tdj,). The Onsager transport equations are 
modified accordingly. The flux J n and \I/ are not of the same tensorial nature (in a first 
approach): the first is a vector defined on the sphere, and the second is a scalar. According 
to the Curie principle, Onsager transport coefficients that couple the two processes should 
not exist. The modification due to the contribution of the electronic relaxation is taken into 
account by a third relaxation term (we still assume T constant) : 

f e 9/2 19/2 

j* _ r dfi 19/2 (g4) 

* = LAfi 

where A/j, is defined as the chemical affinity of the corresponding reaction, or equivalently, 
defined as the pumping force associated to the flux \1/ of spins transferred between the spin- 
polarized electric system and the ferromagnet. In other terms, A/i is the force responsible for 
irreversible spin-transfer. Integrating over the whole ferromagnetic layer, the Fokker-Planck 
equation, Eq. ([55]). rewrites: 

d ~^T = ^ + Ia Lsd AKz)dz = ° (85) 

where the first term in the right hand side ^f- leads to the standard rotational Fokker- 
Planck equation (defined by the second equation in Eqs. (J62j) ) . and the second term is the 
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contribution coming from the electronic relaxation. This equation is the main result that 
defines the irreversible spin-transfer effect for an open system. Eq. ( l83j) will not be solved 
here. The aim of the following developments is to define negative damping and effective 
temperature. 

The dynamic equation is obtained in the same way as in the previous section, by writing 
the flux of representative points on the unit sphere. A term du ^° (not explicit here) related 
to the spin-polarized electronic contribution is added to the Eq. ()63|) and leads to the 
following generalized LLG equation: 



du 
~dt 



du 



•dec 



dt 



g'ux [VV ma9 ) + h'ux 



ux W 



rmag 



(86) 



Since the effe ct of the environment conserves the modulus of the magnetization, i.e. 

telec 



u.du/dt = |116| . the contribution of reduces to the two damping factors, parallel to 
u x VK™ 9 and parallel to tlx (u x \7V ma9 ), and a stochastic force f(t). This necessarily 
leads to the introduction of two parameters a% and a 2 such that : 



^ = -(g' + ai )ux (Vr""'") (// + ,.,) ,7 • 



f(t) 



(87) 



where the coefficients a± and cti can be thought of as negative damping or positive damp- 
ing, depending wether the spin transfer \j/ s ^, integrated over the whole layer with the two 
junctions, is transferred from the electric system to the magnetic system (negative damping) 
or inversely (positive damping). In other words, it depends on the balance of spin accumu- 
lation at the two interfaces of the ferromagnetic layer. Note that as far as the damping 
coefficient s are not explicitly defined, the validity of the argument used above (the Callen's 



argument 



116|) is not restricted to the relaxation and spin accumulation mechanisms de- 



scribed in the first sections of this work, but is much more general. In particular, the equation 
is formally simi l ar to tha t described in t he framework of the exchange torque or spin torque 
theory [20, 
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|. However, beyond the vectorial argument proposed 
above, the proper derivation of the reduction from Eq. (jHHjl to Eq. (|8*7jl (see e.g. [4l| ) is 
still to be performed. 
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C. Effective thermostat 



The Neel-Brown activation laws describe out-of-equilibrium spin systems (dp 7^ 0), and 
are valid for high enough potential barriers kT <C vV ma9 , or long time scales At/r ^> 
1, where At is the measurement time window, and r is the relaxation time scale that 
describes the coupling to the lattice. If the volume v tends to zero, the energy barrier 
decreases down to a value such that vV mag < kTln(At/T ), and the magnetization is at 
" equilibrium" with the lattice for the measurement time window At. The system is no longer 
metastable but superparamagnetic. The equilibrium imposes the condition dp = 0, and the 
statistical distribution is the Maxwell-Boltzmann distribution n = N t exp (— vV ma9 /kT). 
The magnetization behaves like a paramagnetic spin, with the ferromagnetic order parameter 
M instead of the spin s. In the case of 3d metallic ferromagnetic nanostructures (Co, 
Ni, Fe . . . ) of sizes typically around 10 nm radius at room temperature, the system is 
superparamagnetic for time scales of magnetometric measurements (above 10~ 5 sec). The 
system is nevertheless ferromagnetic and follows the Neel-Brown laws if the measurements 
are performed in a shorter time window At (from micro-seconds to nanoseconds in the case 
of 10nm 3 particles invoked above). 

However, if the relevant time window is shorter than the typical ferromagnetic relaxation 
time scale To (At < r ), the precession now governs the quasi-ballistic dynamics, which is 
qualitatively different (because it is not driven by the fluctuations). The collective modes 
measured are that observed with ferromagnetic resonanc e; e.g . dynamics of thermal spin 



waves are observed in GMR structures 1 96 
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98, 



99J, Il00l | . It is no longer activated, 



whatever the potential barrier and the volume v, and the thermalization process vanishes 
at short time scales ("quasiballistic magnetization reversal" regime 9^|). The temperature 
of the system (if any [l23j) is not necessarily the temperature of the lattice T e ff 7^ T: the 
system is decoupled from the heat bath. A similar situation justified the introduction of the 
conc ept of spin temperature T s in the early 50's with the first spin resonance experiments 
|l24| . Note that if the system is also isolated from other sub-systems at comparable time 
scales, has an upper bound, and if the populations (up and down) can be invert ed (like in 



nuclear spin systems), the spin temperature of the system can even be negative 125, 126] 
(but the spin temperature is usually positive and higher than the lattice temperature 126| ). 
In the situation of interest, with spin-polarized currents in 3d metallic nanostructures, 
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the ferromagnetic order para meter is couple d to the lattice through the polarization of the 



electronic degrees of freedom 



105 



127 



128J. Without being coupled to the magnetization 
(e.g. in the non- magnetic side of a junction), the spin-accumulation sub-system relaxes 
toward equilibrium with the relaxation time r s f of some picoseconds (as described in the first 
section with the two channel model). This relaxation time is shorter than the thermalization 
of the magnetization with the lattice tq (nanoseconds). On the other hand, the coupling 
with the ferromagnetic sub-systems corresponds to a relaxation time r s d (comprised between 
the electronic relaxation time r e and r s j) shorter than r : as a consequence, this relaxation 
"thermalize" the ferromagnetic order parameter with the spin-accumulation sub-system, 
that takes the role of the heat bath in place of the lattice. This picture is that schematized 
in Fig n in the first section. 

In the activation regime, it is possible to assume that the spin-accumulation sub-system 
is a reservoir of energy, and that the ferromagnetic order parameter is thermalized with it 
(see Fig. EJ. Providing that the spin-accumulation su b-sys tem is not thermalized with the 
lattice, the zeroth law of thermodynamics is not valid [l 231 ] . and it is possible to identify it 
as a thermostat at temperature T e ff in equilibrium with the ferromagnetic system. 

The equilibrium condition imposes that 1=0 j^J. The entropy production I of the 
spin-dependent electric sub-system was calculated in Sec. II. With the temperature T P ff 
corresponding to the effective equilibrium temperature, the chemical potential writes 38j] 
Afi eff = A/i e// + kT eff ln(n a /n^). 

T eff .l = (~ dA Q° ff + 2eA/i e// ^ 5 J + A/i e// ^ = (88) 

where the Joule heating contribution —^-Jo has been removed because it does not con- 
tribute to the magnetic system and is coupled to the lattice (the whole analysis should also 
include the Peltier effects: energy can also be extracted from the lattice to the magnetic 
system) . 

The entropy production vanishes for the following sufficient condition: 

A/t e// = (89) 



The condition Eq. (J89|) leads to the expression of the equilibrium temperature T< 



*ff- 
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where the inversion of population implies that 5n < 0. This equation is simply the equilib- 
rium Curie- Weiss law that accounts for the paramagnetic behavior of the spin-accumulation 
gUs^n i.e. the first order approximation of the averaging over the Boltzmann distribution 
at temperature T e ff ( 5n is the s-d spin accumulation that would be measured with a lattice 
temperature T e ff). The evaluation of 8n would necessitate the non-equilibrium distribution 
at temperature T to be calculated. This task is beyond the scope of the present review. 
However, according to the evaluation performed below, an energy kT e ff w leV can be ex- 
pected by calculating the ferromagnetic energy under a current of 1 mA due to the spin 
transfer in the internal field of IT of the ferromagnet. 

A fundamental consequence of the existence of the effective temperature is that the 
solution of the stochastic equation of the magnetization is known, and is given by the 
standard activation equation Eq. (|74|) with the effective temperature T e ff instead of the 
lattice temperature T: 

dN 1 9V „ I e ff 

1 NT = -11 Ql) 

09 kT eff 09 2k' eff Trsin(9) 1 ; 

where I e ff is calculated with the Boltzmann distribution with the ferromagnetic energy 

exp(V(6)/kT e ff). Assuming an approximatively constant effective temperature T e ff{9) ~ 

T e ff, the equation is formally identical to Eq. ()74|) so that the Neel-Brown activation formula 

is recovered with T e ff instead of T: 



t{Vi) = r 01 e 
r{9 2 ) = r 02 e 



AV, 



(92) 



AV2 
kT eff 



where Tm an d Tn ? contain the k ' ^ dependence. This behavior is experimentally observed 
|2jJ. Fig. ITU1 shows the typical Neel-Brown activation observed with the sample 
shown in Fig. El(b) and (c). The fit of the mean relaxation time as a function of the applied 
field and the current amplitude injected in the device is performed with the Neel-Brown 
formula with the effective barrier height as fitting parameter (Fig. [TJ3 (a) and (c)). The 
effective barrier height as a function of the current is presented in Fig. E3( c ) an d (d). 

Assuming that the mechanism responsible for T e ff is the spin-accumulation occurring 
at the interface composed of antialigned ferromagnets, the voltage drop due to spin- 
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FIG. 10: Observation of the Neel-Brown activation due to current injection in a Co/Cu/Co 
trilayers (a),(b) and in a Ni nanowire (c),(d). (a) Ratio of the two relaxation times (TLF) and (c) 
relaxation time as a function of the applied field for different currents fitted with the Neel-Brown 
formula and effective barrier height as fitting parameter, (b) and (d) Variation of the effective 
barrier height in Kelvin. In the case of the two level fluctuations, the barrier is measured for the 
symmetric relaxation times (b) (i.e. for different applied fields). The effect of the fluctuations 



: 3 J-e. 



Wegrowe 



are sketched in the insets (dashed lines). Reprint with permission from Ref. 
Phys. Rev. B 68, 214414 (2003) and Ph. Guittienne et al., IEEE Trans. Mag-37, 2126 (2001), 
Copyright@American Physical Society. 

accumulation is approximatively equal to A/z (see Eq. (|18|l): 

AR aa I AR sa I 
kT eff pa — — ; — - oc 93 

/; ln(n a /n 7 ) 2 v ; 

with — Zn(n a /n 7 ) < 1 The proportionality betwe en kT P ff and R sa was observed by dif- 
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131] (with DC measurements, the 



ferent groups in recent experimental investigations 
parameter is the "critical current" I c oc as shown in Fig. El (d)). 

How to estimate the magnetic energy of the spin-accumulation system in usual experi- 
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mental situations, where a current density of some few mA is injected in the nanostructure 



[J l2l Li [s|. I2JJ, |22j, |23j, |2j, |25|? This current corresponds to some 10 16 spin per seconds 
flowing through the interface. If one assumes that 20 % (polarization of the current) of 
the spins are maintained out-of-equilibrium within a typical relaxation time r s f of 10 _11 sec, 
we are left with about 2.10 4 spins that define the magnetization of the spin-accumulation 
sub-system in the volume defined by the corresponding diffusion length. An effect of the 
electric spin relaxation on the ferromagnetic order parameter should consequently be ex- 
pected for a nanostructured ferromagnetic system that is only ten to hundred times larger. 
In an internal field of H int = IT, this energy E = 10 4 fiBH int is of the order 1 eV (beyond the 
Curie tem per ature) in accordance with activation experiments performed on various systems 
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Without current injection, the magnetic order parameter is at room 



temperature and consequently, the hot sub-system is the spin-accumulation system. 



V. CONCLUSION 



An unified thermokinetics approach of both spin-dependent charge carriers and ferro- 
magnetic Brownian motion has been presented in the context of open systems. The spin- 
dependent electronic relaxation is then introduced as a source term in the conservation equa- 
tion of the magnetic moment. This leads to the description of the effect of spin-injection 
induced magnetization switching, or irreversible spin-transfer in an open ferromagnetic layer. 

The description of the spin-accumulated charge carriers is based on the two conduction 
channel approximation, generalized to both intra -and inter- band relaxation. The applica- 
tion of the first and second laws of thermodynamics, together with the conservation laws, 
lead to the spin-dependent transport equations. The relevant Onsager transport coefficients 
are introduced and related to the typical electronic relaxation times. The effect of charge 
conservation and screening is also taken into account. 

On the other hand, the ferromagnetic order parameter is described on an equal footing 
by introducing the conservation laws and the relevant chemical potential, with deterministic 
terms accounting for the effective field, and dissipative terms accounting for the coupling 
to a relevant heat bath. The corresponding Onsager transport coefficients are defined with 
the second law of thermodynamics, and refined with the help of the Onsager reciprocal 
relations. The rotational Fokker-Planck equation, and the Landau-Lifshitz-Gilbert (LLG) 
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equation are then derived within the thermokinetic theory. The Onsager coefficients are 
related to the typical time scales of the ferromagnetic relaxation (r ). In the activation 
regime, the Neel-Brown activation law is deduced. 

In the framework of this description, the generalization of both the Fokker-Planck equa- 
tion and the LLG equation with adding the contribution of spin-accumulation is straight- 
forward in terms of flux of representative points in the magnetization sphere. The negative 
damping appears naturally in order to describe the exchange of spins from the electric sub- 
system to the magnetic sub-system, described as a coupling to an environment. 

Furthermore, the discussion about the different relaxation times shows that the spin- 
polarized current is not thermalized to the lattice in the stationary regime, but is thermalized 
with the spin-accumulation sub-system. The argument is that on one hand the relaxation 
toward equilibrium of the spin- accumulation system (described by r s f. some tens to hundreds 
of picoseconds) is shorter than the thermalization of the ferromagnetic system (described 
by r : nanoseconds). And on the other hand, the coupling between the ferromagnetic order 
parameter and the spin- accumulation sub-system [j s d) is shorter or equal to the r s /. 

An effective temperature is then derived in the activation regime through the entropy 
production, and leads to the derivation of an effective Neel-Brown relaxation process due to 
current injection, that is experimentally observed. 
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APPENDIX A: MICROSCOPIC APPROACH AND THERMOKINETIC COEF- 
FICIENTS 



1. Relation between L and the electronic relaxation times 

Let us consider a simple interface between two metals. Far from the interface, the Ohm's 
law is recovered : the chemical potentials of the channels are identical and the electric 
distribution is that of equilibrium n aQ and n 70 - In the following, we assume that the charge 
transfer between the two channels can be described by the following relation: 

fAn a (x) + gAn^x) = (A92) 
The case / = g = 1 describes the local electrical neutrality. We have: 



An a (x) = (/i cM - fj?)N a (E F ) 

An 7 (x) = (/i^ - /jP)Nry(E F ) 



(A93) 



7 



where /i is the chemical potential in the absence of charge transfer and Hch,a and Hch. 
are the purely chemical potentials of the channels (without transfer fi c h,a = Hch,-y = A*q) ; 
N aa (Ep) is the density of states at the Fermi level. The separation between the electric 
potential and the purely chemical potential writes: 



(A94) 



Ha = Hch,a + eV 

Hi = Hch,-y + eV 

where V is the local electric potential. 
Relation (|A"92j) gives 

{Hch,a ~ H°)fN a (E F ) + {fx ch>1 - fj,°)gN 7 (E F ) = (A95) 

Note that the expressions of An a 7 comes from the fact that at zero Kelvin: 




N a {E)f{E)dE= r""* N a {E)f{E)dE 



OO 



(A96) 



N 7 (E)f(E)dE = / N 7 (E)f(E)dE 
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From these relations we deduce: 
An a = N a (E F ) S/i chja 



gN a (E F )N 1 (E F ) 



fN a (E F )+gN^E F ) 
9 I 7 fN a (E F )+gN 7 (E F s 



(A97) 



ch,a 



where 5// a7 = /i ch)£n , - // - 

Introducing the global transfer rate T Q ^ 7 (resp. T 7 _, a ) of the channel a to 7 (resp. 7 to 
a), the charge conservation between the two channels writes (e < 0): 



dn a (x) ldJ a (x) 



dt 

dn y (x) 



dx 



T a ^ y (n a (x) , n 7 (x) ) + T^ a (n a (x) , n y (x) ) 



ldJ y {x) 



a— >7 



n a (x) , n 7 (x) ) - T 7 ^ a (n Q (x) , n 7 (x) ) 



<9t e <9x 

that leads, in the stationary regime ^.t-O) _ q ? ^ foe following relations: 



(A98) 



dJgjx) 

dx 
dJj(x) 

dx 



= -eT a _> 7 (n a (x),n y (x)) + eT^ a (n a (x), n 7 (x)) 
= +eT a ^ y (n a (x), n 7 (x)) — eTj^ a (n a (x), n y (x)) 



(A99) 



The Taylor expansion to the leading order of the transfer rates T Q ,_ +7 and T 7 _> a , around 



equilibrium gives: 

r oj a (x) 

dx 



dT _> dT _> 



dn 
dT 

V ± 7 — >Q 



dn^ 
dT^ 



+ eT 7 _«, n 7 ) + e^An a + e^A« 7 



dJ 7 (x) dJ a (x) 



dx dx 

At equilibrium, the current of each channel is conserved, so that: 

-eT a ^ 7 «,n 7 ) + eT 7 ^ Q «,n 7 ) = 
Defining the electronic relaxation times: 



(A100) 



(A101) 



1 ^ (9(T 1 Q ,_ >7 Tly^Q,) 

1 d(T^—> a T 1 Q ,_ >7 ) 



(A102) 
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we have 



dJ a (x) 

dx 
dJ^(x) 



—e 



dJgjx) 

dx 
dJ 7 (x) 

dx 



N a (E F )N^(E F ) 



—e- 



+ e- 



Aa, 



dJ a (x) 



dx 



fN a (E F )+gN 1 (E F ) V 
dJJx) 



9 f w 



dx dx 

The above equations can be rewritten in the following form : 



(A103) 



(A104) 



dx 

dJJx) _ . 



c/1,7) — ~\~L ~ P"y) 



(A105) 



where the coefficient Onsager transport coefficient L is related to the electronic relaxation 
times by the following relation: 



L = e 



N a (E F )N^{E F ) 
fN a (E F )+gN^(E F ) V 



9 + / 



The form Eq. (|A105|) is that of Eq. © deduced from the thermokinetic approach in Sec. 
I, where the coefficient L is the Onsager coefficient defined in the third equation of Eqs. 

dnj. 



2. Determination of / and g 



We have, for the stationary regime: 



d(J a (x) — Jy(x)) 
dx 



-2L(fi a - /i 7 ) 



(A106) 



Furthermore, the local Ohm's law applied to each channel leads to the following equations: 



J a 



a a d\i a (x) 

e dx 
<r 7 dfj, 7 (x) 

e dx 



dJ a (x) 

dx 
dJ^(x) 

dx 



a a d 2 fi a (x) 

e dx 2 
<r 7 <9 2 /i 7 (x) 

e dx 2 



(A107) 
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where we assume that the conductivities are constant. 
From (IA105D and (IA107I) . we get: 
d 2 ii a (x) eL 



9 /x 7 (x) _ eL dx 2 \cr Q a 7 

that leads to the well-known diffusion equation of the chemical potential, that describes the 
spin- accumulation process: 



from which the spin-diffusion length l s f is deduced: 



1 , / 1 1 



From Eq. ()A108|) we have the differential equations: 



p -e^l— + -] ( A11 °) 



d 2 ^ a (x) eL ^ o-jd 2 Afi 

dx 2 <y a <y t dx 2 

<9 2 /x 7 (x) eL a a d 2 A/j 



(Alll) 



dx 2 o~j ^ o t dx 2 

3. Charge distribution and screening 

Separating the electric contribution from the chemical contribution, the electrochemical 
potential writes : 

/Vr = Vch,<xa + eV (A112) 

so that: 

dx 2 ~ dx~ 2 + 6 dx 2 ^ I dx 2 o t dx 2 6 e 

d 2 fj 11 (x) d 2 fi cht y(x) d 2 V(x) I d 2 fi chr/ (x) a a d 2 A^i 2 An Q + An 7 



dx 2 dx 2 dx 2 K dx 2 a t dx 2 e 

(A113) 

where the Poisson equation has been introduced 

d 2 V(x) A^ + An, (An4) 



dx 2 
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The equations rewrite: 



f d 2 8fi ch ^(x) _ 2 N a 5y,ch, a + NjSf^chsy = ^ 

e a t I 2 



dx 2 



sf 



e a t I 2 



These relations with 



lead to 



dx 2 



fl a fij fl>ch,a l^chj-y d^ch,a ^ch,-) 



sf 



(A115) 



(A116) 



' d 2 Afi cKa {x) A/i chjC 



dx 2 I 2 
d 2 Afi chn (x) Afi chn 



= Afj, 

= Afl 



2^7 



sf 



dx 2 I 2 
where we have introduced the screening length : 

1 2 A a + A 7 



a t l 2 



+ e 



2 N a 



sf 



(A117) 



p - e (AIM) 

The solution of the equation is composed by a solution of the equation with zero right 
hand side (homogeneous solution) and a particular solution. 



1. Solution for Afi^ a 
Homogeneous solution 

A^ = Aexpfi)+Bexp(-f) 
Particular solution A//^ r * = p a A/x 



(A119) 



Afi Afj, ( a y 2 A 7 



72 72 



"sf 



Pa = 



1 1 

T 5 ¥ 



(A120) 
(A121) 
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2. Solution for A/i c ^ 7 
Homogeneous solution 

A/i™ 9 = A'exp (|) + 5'e*p (-y) (A122) 
Particular solution: A/i^^ = p 7 AyU 

A^ A// /-<T Q 2 A" a \ /A10Q^ 

^p 7 _-_=^— + e -)A M (A123) 

gg _|_ e 2 Af a 

=► P7 = (A124) 

The general solutions satisfying the condition 

Aii ch>a - Afx ch ^ = A/i (A125) 

correspond to 

A = A' B = B' (A126) 
Inserted in the expression of the charge conservation 

fAn a + (?An 7 = (A127) 

we have 

fN a (Aexp (y) + Bexp (-y) + p a A/i) + #7V 7 (Aexp (y) + Bexp (-y) + p 7 A^) = 

(A128) 

(/iV a + gN 7 ) {Aexp (y) + Bexp (~)) + (fN a p a + gN^)Afi = 0, for all x 

(A129) 

f'N a p a + #A" 7 p 7 = and A = B = (A130) 



A solution writes: 
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f = N, 



Qg _ e 2 N a 
a t f sf e 



(A132) 



4. function of l a ,l^ and / 



The relation between L and the electronic relaxation times has been found to be: 



= N^WEr) tj_ _J_\ (A133) 
Inserting the expression of / and g obtained in the previous paragraph, L becomes: 

1 / (TaN-y . o-,N a \ 2N a N~, ( l 



L = e — j j (A134) 



FT ' 2 

s/ 



Furthermore, according to Eq. ([AllOjl . the coefficient L can also be written in the 
following form: 



L = (A135) 



ea t l sf 

Both results lead to the equation: 



cf a N^ 
&tl 8 * y ^ — >ck 



+ 



^2NaNj_ ( 1 



— »7 



+ 



"T^ — >a 



a t l 2 



sf 



Let us define the typical diffusion length per channel l a et L such that 



e 2 AL 



Eq. (|A136|) rewrites: 



1 a a 


r 0-^7 




t e 2 N I 2 


Ty—Hx 


e 2 l 2 


I rr- 


a^N a 


(J-y(J a 






e 2 l 2 



(A136) 



(A137) 



e 5 ^ 



IT,, ^- / J 1 1_ 

^7 



N a ll 



+ 



AT 7 i 7 



s/ 



<T a (T'Y 



+ 



aJ2 I sf 



1 1 

l 2 + I 2 



(A138) 
(A139) 
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a. Limits: In metals, the screening length is much smaller than the diffusion length of 
both channels: 



1. l/lcty < 

0t 



Ny N a , ., 



N a + N y \<Jyll acPJ Sf 



The other limit gives: 



Ny/{N a + Ny) 1 



J 2 



N a /(N a + jV 7 ) 1 
(7 7 /(cr a + (7 7 ) /2 1 er a /(er a + er 7 ) Z 7 



+ 



(A140) 



2. Z/Z 



CV7 



OO 



I 2 I 2 I Sf 



1 



s 



(A141) 



But this second limit is not expected in usual materials. 



APPENDIX B: ONSAGER MATRIX 



The aim of this Appendix is to derive the Onsager matrix ([28)1 . (and as a particular 
case (jllj) ) on the basis of the first and second laws of thermodynamics. In a typical one 
dimensional junction the layer is decomposed into Q parts, defining the sub-system T, k , 
which is in contact to the "reservoirs" E fc_1 and S fc+1 . The sub-systems E fc , is then an open 
system which exchanges heat and chemical species with its left and right vicinity layers. 
Furthermore, the populations (NK, A 7 ^, N$+) and spin down (N^) are not conserved due 
to transitions from one channel to the other. 

In this picture, the states of the sub-system E fc are described by the variables 

(S\N^N^N k dp N k dL ) (B142) 

where S is the entropy. The internal variables ^ and ty s d must however be introduced 
in order to take into account the relaxation processes occurring respectively between the 
two s-spin channels, the two d-spin channels, and the s-d relaxation . 

Let us define the heat and chemical power by P$ (the mechanical power is zero as long 
as the action of the magnetic field on the charge carriers is neglected). The first law of the 
thermodynamics applied to the layer S fc gives 



dE k 
~dT 



P. 



P 



k^k+1 



(B143) 
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Introducing the canonical definitions T k = ||£ and ji k s± = J^r, ^ d ± = -§^~ the energy 
variation is: 



dE k _ Tk dS k k dN k , dN k , dN k , dN k , 

For the sake of simplicity, we limit our analysis to the isothermal case, T k = T. The 
entropy variation of the sub-layer is deduced from the two last equations, after introducing 
the conservation laws and after defining the polarized currents 51 \ = I s ± — 1^, 51 \ = 
I si - I dl , and the currents I L = I sl + I dl , I s = I s] + I sl , 



T^- = P*" 1 "** _ p k^k + i _ 1 AfJ k (sjk-i^k _ 6I k^k+i + ^ _ 2 ^ 

\» k (it 1 ^ - i^ k+1 - - \ M (sit' 1 "" ~ 5I V k+1 - 2 ^ - 



"2^ 
1 



2 -/x{ (/J- 1 -* - + #*) (B145) 

where we have introduce the chemical potentials /xj = fi k ^ + A*^/2, //J = + /x^, 

and £/ie chemical affinities of the reactions, defined by A/iJ? = //^ — [JL k j2 = — §fr, AyuJ = 
,.k In ,.fc _ _ d£* 

set 

The entropy being an extensive variable, the total entropy variation of the system is 
obtained by summation over the layers 1 to Jl where the layer 1 is in contact to the left 
reservoir R l and the layer Vt is in contact to the right reservoir R r . 

The total entropy variation is: 

T— = [P] Rl ^ 1 - [p} Q ^ Rr 
dt 

n n 



+ E \ - SI s-^ k + E \tfr x - A) Ik ^ k 

k=2 k=2 

+ E^(^r-M) ^ k + 

k=2 

si 1 n n 

-M~ l - »\)i k ^ k + E A ^ (Bi46) 



2 

k=2 k=l k=l 



where [P]^^ 1 and [p]^-*^ i n the right hand side of the equality stand for the heat and 
chemical transfer from the reservoirs to the system E. 
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The entropy variation takes the form 

= ^FiX* + P ext (t) (B147) 

i 

where Fi are generalized forces and X 1 are the conjugate generalized fluxes. The variation of 
entropy is composed of an external entropy variation P ext (t)/T and by an internal entropy 
variation I. 

By applying the second law of thermodynamics I> we are introducing the kinetic coef- 
ficients Lij such that X = £V Fi \^2j LijF^j . By identification with the expression (|B146|) . 
the kinetic equations are obtained, after performing the continuous limit, 



Jos 




L S s 


















-8 lis 

dz 


Joi 




L[ S 


L ll 
















— 8AfJ, s 

dz 


5Jf 










LgsSs 












—0/i(5 
dz 


5Jf 










Td 


Td 










dz 

























A/i s 





















ji 
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The kinetic coefficients are state functions; 



matrix is positive : 
laxation processes (j 



Lij(S k , N+, N*) and the symmetrized 
Lji + Lij}^y > 0. The coefficients L int refer to the internal re- 
According to Onsager reciprocity relations, the kinetic coeffi- 



cients are symmetric or antisymmetric = ±Ljj. The coefficients are known from the 
two-channel model for the conductivity. The two last equations concern the internal (L int ) 
"density" variables ip s and ip s d defined by \l/ fc = ip(z)dz. Due to the Curie principle, there 
is no coupling between spin polarized transport processes and the electronic transitions (the 
scalar process is not coupled to vectorial processes). 
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